f=@(x,y)exp(-x.^2/2).*sin(x.^2+y);
J=integral2(f,-2,2,-1,1,'RelTol',1e-20)